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Abstract 

The shear dynamics in Bianchi I cosmological model on the brane with 
a perfect fluid (the equation of state is p = (7— is studied. It is shown 
that for 1 < 7 < 2 the shear parameter has maximum at some moment 
during a transition period from nonstandard to standard cosmology. An 
exact formula for the matter density /1 in the epoch of maximum shear 
parameter as a function of the equation of state is obtained. 

In recent several years a lot of effort has been done on the idea that our 
Universe is a boundary of a space-time manifold with a larger number of di- 
mensions jjj. The cosmological evolution of a three-brane in five dimensional 
space-time become the matter of intense investigations. The corrections to Ein- 
stein equations on the brane arising from the influence of a bulk geometry on 
the brane were obtained The corresponding cosmological dynamics was 
studied in detail for a FRW Universe. It was shown that the dynamics of the 
early brane Universe can differ significantly from the standard scenario ||. On 
the other hand, the standard FRW cosmology for the late Universe (when the 
matter density on the brane /i become negligible in comparison with the brane 
tension A) was recovered 

Later on this analisys has been done for an anisotropic Bianchi I metric on 
the brane. The matter source was chosen in the form of a massive scalar field 
H or perfect fluid with the equation of state p = (7 — l)/x ](| (see also [Q 
where some exact solutions were found). One of the most interesting results 
obtained in this way is that the initial singularity, being anisotropic in ordinary 
cosmology (with the exception of the maximally stiff fluid case ||), becomes 
isotropic in the brane cosmology for 7 S (1, 2)|| ^j. Since during the evolution 
of the brane Universe its dynamics approaches the standard one, later stages 
of the Universe are isotropic for 7 < 2. So, the intermediate stage of a rela- 
tively high anisotropy can exist if the matter filling the Universe has a positive 
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pressure. Intuitively we could expect that the anisotropy maximum is reached 
sometimes in the transitional period when the matter density is comparable 
with the brane tension. In this paper we calculate the ratio fi/X during the 
maximum anisotropy epoch. We will see that for the radiation-dominated Uni- 
verse the above suggestion appears to be true, though for an arbitrary 7 from 
the interval (1,2) (we study only this type of matter in this paper) the epoch 
of maximal anisotropy can differ significantly from the period when \x ~ A. To 
conclude this introductory part we remind that current experimental limit on 
the brane tension is A > (lOOGeV) 4 §. 
The field equations on the brane are M 
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Here h is the 5-dimensional gravitational constant, k is the effective 4-dimensional 
gravitation constant on the brane. Below we assume that the effective cosmolog- 
ical constant on the brane A = 0. Bulk corrections to the Einstein equations on 
the brane are of two forms: there are quadratic energy- momentum corrections 
via the tensor S^ v and nonlocal effects from the free gravitational field in the 
bulk. The energy-momentum corrections are given by 
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where T = T M M . If we define w M as the 4- velocity comoving with matter, the 
nonlocal term has the form 
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where = + u tl u u . Here 

U = --k 2 \£ 
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is an effective nonlocal energy on the brane, 
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is an effective nonlocal anisotropic stress and 

is an effective nonlocal energy flux. 

There are also conservation equations for the brane energy-momentum tensor 
SJ^T^v = 0. Full description of dynamical and conservation equations see in Q. 

We will study Bianchi I model on the brane with the metric 



ds 2 
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The conservation equations in this case reduce to 



fi + e^+p) = o, (5) 
u + \eu + ct^tv = °, ( 6 ) 

£>"7V = °i ( 7 ) 

where is the volume expansion rate, and <Jn V is the shear. The nonlocal energy 
flux Q vanishes in this case identically and we have no evolution equation for 
T\ v which is a bulk degree of freedom and can not be predicted from the brane 
[U |). In terms of the mean scale factor a — (aia^aa) 1 / 3 the expansion rate 
9 = 3H = 3a/ a. 

The Raychaudhuri equation on the brane is H 



6 + \e 2 + v^a^ + ^ 2 (p + 3p) = -\k 2 {2h + 3p)^ %L , s) 
and Gauss-Codazzi equations are (f| 



0>„ + 0O>„ = —^rV^u, (9) 
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Due to the presence of the nonlocal stress V^v this system of equations is not 
closed. In previous papers || [| this difficulty was avoided by choosing U = 0. 
This condition, however, is sufficient, but not necessary for closing the system. 
The full situation can be described as follows. 

The most restrictive condition is V^v = 0. In this case it is possible to solve 
the equation for the shear separately for components of the shear tensor. The 
less restrictive condition a^V^u — allows us to find the dynamics of shear 
scalar 'er M „ leaving components of the shear tensor undefined. 

Taking into account Eq. (6), we can see that both these assumptions lead 
to Friedmann-like behavior of the " dark radiation" term 

with C = Const. 

The condition U — automatically leads to u^V^v — (see Eq. (6)). How- 
ever, the latter condition is more general. We assume further that <r A ""P AI „ = 
is satisfied. In the opposite case the nonlocal anisotropic stress changes signif- 
icantly the dynamical equations for shear and "dark radiation". Though even 
in this case there are some possibilities for closing the system (8)-(10) Q the 
resulting dynamics can be very different from studied in this paper and should 
be specially investigated. 
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It should be noted that though there are no restriction to choose V^u on the 
brane, we do not know whether a particular chose is consistent in a bulk. The 
general form of bulk anisotropic metric is not known yet (the first approach to 
this problem see in [|ll|). We leave the problem to check that V^ v chosen on the 
brane is consistent with the full 5-dimensional metric to future investigations. 

It is convenient to introduce the deceleration parameter q as 

da 
a 1 

or, equivalently, 

H = -(l + q)H 2 . (11) 
This equation in combination with Eq.(8), which can be written in the form 

9 2 9 k 2 ll , k 2 ll 2 , 2U , , 

H = ~ H ~ r ~ - 2 > - -£r& -»-*x> (12) 

gives the expression for q. To obtain this expression in more convenient form 
we introduce dimensionless variables 

s2 = i- < 13 > 

Substituting Eqs.(13)-(16) into (12) we obtain 

(7 = 2S 2 + i(37-2)^ + (3 7 -l)fl A + r! w (17) 
and the constraint equation (10) in the form 

i = s 2 + n A + n M + n w . (is) 

With a new time variable r defined as 4^ = -h the dynamical equation for 
the shear takes the form as in standard Bianchi I cosmology j|] 

(£ 2 )' = 2( (? -2)£ 2 (19) 

where ' is the derivative with respect to r. We remind, however, that unlike 
standard cosmology we can not determine shear components separately unless 
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Vp.u = 0. The Eqs. (17)-(19) take the form obtained in in the particular case 
U = 0. Using (19) we can easily find that the extremum of the shear parameter 
E corresponds to deceleration parameter q = 2. Expressing E from (18) and 
substituting into (17) we find that 

g = 2 + |n jU ( 7 -2)+30 A ( 7 -l)-n w . (20) 

This leads to the following equation satisfied at the moment when E reaches its 
extremum: 



n w = -n M (7-2) + 3ft A ( 7 -i). (2i) 
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First let us assume U = 0. In this case it is possible to obtaine the result in 
a very compact form. The equation (21) yields now 

fiA(7-i) = ~n„(7-2) 

Since Q \ falls more rapidly than fl^ , this extremum of shear is indeed its max- 
imum. 

Remembering the definitions of fl\ and we can write that = 2A/ /i. 

Using this formula we finally get the ratio of matter density and the brane 
tension at the moment of maximal shear parameter for a given equation of 
state: 

For the most interesting case of radiation-dominated matter (7 = 4/3) the 
maximum of E corresponds to fj, = 2A. 

In the case of the maximally stiff matter (7 = 2) the isotropisation does 
not take place even in the standard cosmology, so for 7 close to 2 the shear 
parameter can reach its maximum at an arbitrary late time with an arbitrary 
small ratio /i/A. In the opposite case of a very low pressure the isotropisation 
can begin at a very high matter density. If the pressure is exactly equal to 
zero (7 = 1) or it is negative (7 < 1), the shear parameter decreases from the 
beginning. In this case there is no maximum of E and formula (22) become 
inapplicable. 

If, as often done, we write the equation of state through the parameter w as 
p = w/i, Eq. (22) takes even more simple and transparent form 

/i 1 — w 



A w 

w G (0,1). 

In the case U 7^ it is necessary to use the full form of (21). It can be 
rewritten as a quadratic equation for /i/A 

£b-l) + ^- 2 > = ^ (23) 
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but with righthandside decreasing with time. For positive hi (otherwise the 
brane Universe can recollape without reaching the isotropic stage ]Io|]) Eq. (23) 
has only one positive root, corresponding to maximum of the shear parameter. 

As hi decreases like radiation (hi = C/a 4 ), the result obtained for radiation 
matter without hi can be easily generalised. For radiation-dominated Universe 
it is possible to write 

U 

with some dimensionless constant parameter a. Eq.(23) gives now the matter 
density at the time of shear maximum as [i = A(2 + 12a). In general, for positive 
hi the isotropisation time shifts to an earlier time than for hi = 0. In particular, 
in the case of the maximally stiff fluid the shear parameter has the maximum 
at 

Clearly, /i — > with hi — > 0, reflecting the fact that there is no isotropisation for 
the maximally stiff fluid without "dark radiation". 

We investigated the shear dynamics in Bianchi I cosmology on the brane. 
Unlike ordinary General Relativity scenario (where the shear parameter mono- 
tonically decreases from Kasner initial singularity to isotropic future attractor 
unless 7 = 2) this dynamics for 7 € (1,2) describes the evolution of the Universe 
from isotropic singularity to isotropic future attractor through an intermediate 
anisotropic stage. It was shown that assuming a^V^v = on the brane it is 
possible to find the matter density in the epoch of maximum shear parameter 
analytically. A more restrictive condition V^v = on the brane allows us to 
separate the evolution of shear tensor components. An important problem for 
future development is to find 5-dimensional bulk metric for a Bianchi I brane- 
world and check the consistency of these forms in the bulk. 
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